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A rigorous microscopic theory for the description of quantum-transport phenomena in systems 
with open boundaries is proposed. We shall show that the application of the conventional Wigner- 
function formalism to this problem leads to unphysical results, like injection of coherent electronic 
states from the contacts. To overcome such basic limitation, we propose a generalization of the 
standard Wigner-function formulation, able to properly describe the incoherent nature of carrier 
CN| ■ injection at the device spatial boundaries as well as the interplay between phase coherence and energy 

relaxation/dephasing within the device active region. The proposed theoretical scheme constitutes 
' the first rigorous derivation of the phenomenological injection model commonly employed in the 

simulation of open quantum devices. 

CLi' PACS numbers: 72.10.Bg, 85.30.-z, 73.40.-c 
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Present-day technology pushes device dimensions toward limits where the traditional semiclassical transport the- 
ory can no longer be applied, and more rigorous quantum-transport approaches are required However, in spite 
of the quantum-mechanical nature of carrier dynamics in the core region of typical nanostructured devices — like 
semiconductor superlattices and double-barrier structures — the overall behaviour of such quantum systems is often 
f-^. ■ the result of a complex interplay between phase coherence and energy relaxation/dephasing ||, the latter being pri- 
marily due to the presence of spatial boundaries Q . It follows that a proper treatment of such novel nanoscale devices 
requires a theoretical modelling able to properly account for both coherent and incoherent — i.e., phase-breaking — 
processes on the same footing. To this end, a generalization to open systems — i.e., systems with open boundaries — 
of the well known Semiconductor Bloch Equations (SBE) || has been recently proposed ||. However, the theoretical 
analysis presented in j(| is primarily related to the interplay between phase coherence and energy relaxation within 
the device active region, and — apart from its abstract formulation — no detailed investigation of the carrier-injection 
process (from the electrical contacts into the device active region) has been performed so far. 

Aim of the present Letter is to provide a rigorous quantum-mechanical description of the coupling dynamics between 
the device active region and external charge reservoirs, able to account for the semi-phenomenological injection models 
commonly employed in state-of-the-art simulations of realistic one- and two-dimensional open quantum devices [M. 
Among such simulation strategies it is worth mentioning the approach recently proposed by Fischetti and co-workers 
By denoting with f a the carrier distribution over the electronic states a of the device, the transport equation 
^ \ proposed in || is of the form: 

ifa = y(W aa ,f a ,-W a , a f a ) + ^- A . (1) 



c3 ; 

Here, denotes the equilibrium carrier distribution in the contacts while r Q can be regarded as the device transit 
time for an electron in state a. As anticipated, in spite of a rigorous treatment of the scattering dynamics, the last 
(relaxation-time-like) term describes carrier injection/loss on a partially phenomenological level and does not depend 
on the real position of the device spatial boundaries. 

In order to provide a fully microscopic real-space description of the carrier-injection process, we shall start rivisiting 
the theoretical approach proposed in |(| . By applying to the conventional SBE for a closed system || |l^, |ll| 

d Paia2 = J2 L ^ M.P°M ( 2 ) 



dt' 



a, a 



the usual Weyl- Wigner transform 



one gets: 
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,k)= / dv'4> ai \v+ r -\ — T 4>* (r-V) (3) 



j t f(r, k) = J L(r, k; r', k')/(r', k')dr' dk' , (4) 
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FIG. 1: Schematic representation of the device active region sandwiched between its electrical contacts (a) and of the 
corresponding U boundary-condition scheme for a one-dimensional system. The latter implies, in particular, the knowledge of 
the incoming Wigner function f(zb,k), i.e., f(zy e ^,k > 0) and /(^right'^ < 0). 



where 

/( r - k ) = ^2 u aiCt2 (r,k)p aia2 (5) 

is the well-known Wigner function |l| [l2| while the quantity L(r, k; r', k') is the Liouville operator in the new phase- 
space representation r, k. 

For a closed system, the Wigner function / is defined for any value of the real-space coordinate r and its time 
evolution is fully dictated by its initial condition. In contrast, for a system with open boundaries / is defined only 
within a given region f2 of interest and its time evolution is determined by the initial condition at t inside such 
region as well as by the value of its incoming part on the boundary of the domain at any time t' > t . Such 
boundary-condition scheme jl3| — also referred to as "U scheme" — is depicted in Fig. [l]. More specifically, in order 
to properly impose the desired spatial boundary conditions, let us rewrite Eq. (^) as 

j t f(r, k) = J L(r, k; r', k')/(r', k')dr'dk' + S(r, k) (6) 
by adding and subtracting a source term S(r, k) = u(k)/ & (k)<5(r — Tb) [|l4j. The latter is compensated by a renormal- 
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ization AL = L — L of the Liouville operator L in (|4|) |15j : 

AL(r,k;r',k') = -v(k)6(r - r b )S(r - r')5(k - k') . (7) 
By applying the inverse of the Weyl-Wigner transform in (||) to Eq. (^) we get: 

~j^Pa 1 a 2 = ^ t L aia2 a i^ a ^p a i^ a i^ + S ai a 2 ■ (8) 



ct , a 



l"2 



Equation (g) is the desired generalization to open systems of the SBE in (|2|) [g_6| . 

In order to validate the theoretical approach presented so far, we shall focus on a very simple semiconductor 
nanostructure: a single-barrier equidistant from the device contacts (see Fig. ||). As basis states a we adopt the 
scattering states of the device potential profile; moreover, to better identify the role played by carrier injection, we 
shall neglect all other sources of energy relaxation/dephasing in the device active region, like carrier-phonon and 
carrier-carrier scattering. Under these assumptions, Eq. (^) in steady-state conditions reduces to: 

^7 ( e cti £a 2 ) Pa 1 a 2 ~ ^ ] ^-l-'aia 2 ,a l 1 a' 2 Pa' 1 a' 2 — S aiOC2 , (9) 
a{a' 2 

where e a denotes the energy of the scattering state a and AL QlQ2jQ ' iQ 7 is the renormalization term in (|t]) written in 
the density-matrix representation a±a2- 

Figure ^| shows results for the single-barrier potential profile when carriers are primarily injected from left. Here, 
the simulated real-space charge distribution obtained from the phenomenological injection model in Eq. ([l]) (dashed 
curve) is compared to that of the microscopic model in (^|) (solid curves). As we can see, the two models give 
completely different results. The phenomenological model gives basically what we expect: since we have significant 
carrier injection from left only and since the potential barrier is relatively high, the carrier distribution is mainly 
located on the left side. In contrast, the microscopic model gives an almost symmetric charge distribution. In order 
to understand the origin of this unphysical result, let us focus on the nature of the source term in (|^). Contrary to the 
phenomenological injection/loss term in (Q), the latter is intrinsically non-diagonal, i.e., the injection of a carrier with 
well-defined wavevector k is described by a non-diagonal source contribution S aia2 . In other words, we inject into 
the device active region a coherent superposition of states a\ and ct2, in clear contrast with the idea of injection from 
a thermal — i.e., diagonal — charge reservoir. More specifically, in this case the generic scattering state a on the left 
comes out to be an almost equally weighted superposition of +k and — k: 4>k{z) = ake lkz +bke~ tkz . This, in turn, tells 
us that the generic plane-wave state k injected from the left contact is also an almost equally weighted superposition 
of the left and right scattering states. This is why the charge distribution (solid curve in Fig. ||) is almost symmetric: 
any electron injected from left couples to left as well as to right scattering states. The anomaly of the microscopic 
model is even more pronounced if we look at the carrier distribution in momentum space (see inset in Fig. |^) . While 
for the phenomenological model (dashed curve) we get a positive-definite distribution showing, as expected, the two 
symmetric wavevector components of the scattering state, the microscopic result is not positive definite; this tells us 
that our boundary-condition scheme does not ensure that the solution of Eq. @ is a "good" Wigner function. 

The scenario previously discussed is highly non-physical; it can be mainly ascribed to the boundary-condition 
scheme employed so far, which implies injection of plane- wave electrons [see source term in Eq. (||)], regardless of the 
shape of the device potential profile. This is an intrinsic limitation of the conventional Wigncr-function representation 
r, k. It is then clear that, in order to overcome this limitation, what we need is a boundary-condition scheme realizing 
diagonal injection over the scattering states a of the device potential profile. 

To this end, in this Letter we propose a generalization of the Wigner-function formulation considered so far. The 
key idea is to extend the Weyl-Wigner transform in (||) from the k to a generic basis set {|/3}} according to 0: 



ufX (r) = 1 dr'^ ai (r + X ft (r + ^) X 02 (r - ^) C 2 (r - . (10) 
In analogy to (|^) , our generalized Wigner function is given by 01| : 

/ ftft (r) = 5]«;;;;(rK ia2 . (11) 
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FIG. 2: Comparison between the real-space charge distribution obtained from the phenomenological injection model in ([j]) 
[n(r) = f a \4>a( r )\ 2 — dashed curve] and the microscopic model in Q [n(r) = ^ a a p aiQ12 ai (r)0* 2 (r) — solid curve] 
for a GaAs-based single-barrier structure (height V a = 0.5 eV and width a — 4nm) equidistant from the electrical contacts. 
In this room-temperature simulation, due to a misalignment A/i = 0.2 eV of the left and right chemical potential, carriers are 
primarily injected from left. The corresponding charge distribution in momentum space is also reported in the inset (see text). 



By adopting as basis states \(3) the scattering states of the device potential profile and assuming a diagonal source 
term of the form 



(12) 



its equation of motion will be 



^//3i/3 2 ( r ) = / dr '^/3i/32,/3^ 2 (i-,i-')//3j/3 2 (r') + S 01 p 2 (r) (13) 

with a renormalization ALp 1 p 2 pi^{r, r') given by 

- VfaSfa^SfrfawJiT - n)S{r - r') . (14) 

We stress that now the source term in (|l^) describes diagonal injection over the scattering states (with velocity vp), 
as requested. Indeed, if we now integrate Eq. (U3) over the real-space coordinate r, we get again the density-matrix 




equation in (g), but with a diagonal source term S ai a 2 = v ai f^ 6 ai a 2 and a much simpler — i.e., partially diagonal — 

renormalization term A.L aia2iQ ,' iQ y = — v ai 5 ai a 2 Ua 1 ial (rt>). In the scattering-free case, the stationary solution is again 
described by Eq. (||). However, due to the diagonal nature of the new source term as well as of the partially diagonal 
structure of AL, Eq. (|^) has now a diagonal solution: p ai a 2 = J "a/a^j- More specifically, the diagonal density-matrix 
elements f a obey the following steady-state equation: 

T~aa' fa' = f a (15) 

with T aa i = (ri,). Equation ( [l5| ) is semiclassical in nature, i.e., it involves diagonal density-matrix terms only [ p)| . 

Figure [3] shows again results for the single-barrier potential profile previously considered. Here, the simulation 
based on the phenomenological injection model in Eq (Q) (dashed curves) is compared to that of the new microscopic 
model in (|l^) (solid curves). As we can see, the highly non physical behaviours of Fig. || (solid curves) have been 
completely removed. Indeed, the momentum distribution in the inset is always positive-definite and the two models 
exhibit a very similar behaviour. We find relatively small deviations close to the device spatial boundaries, which 
can be ascribed to the interlevel injection coupling T aa i [see Eq. (p"5|)], not present in the phenomenological injection 
model. 

In conclusion, we have proposed a rigorous description of quantum-transport phenomena in systems with open 
boundaries. Our analysis has shown that the conventional Wigner-function formalism leads to unphysical results, 
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like injection of coherent superpositions of states from the device spatial boundaries. This basic limitation has been 
removed by introducing a generalization of the standard Wigner-function formulation, able to properly describe 
the incoherent nature of carrier injection. The proposed theoretical scheme constitutes a rigorous derivation of the 
phcnomcnological injection models commonly employed in the simulation of open quantum devices. 

We are grateful to Aldo Di Carlo, Massimo Fischetti, Rita Iotti, Carlo Jacoboni, Paolo Lugli, and Paolo Zanardi 
for a number of stimulating and fruitful discussions. 
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